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' The nonlinear amplitude modulation of longitudinal dust lattice waves (LDLWs) propagating in a 

I dusty plasma crystal is investigated in a continuum approximation. It is shown that long wavelength 

, LDLWs are modulationally stable, while shorter wavelengths may be unstable. The possibility for 

^SJ ' the formation and propagation of different envelope localized excitations is discussed. It is shown 

that the total grain displacement bears a (weak) constant displacement (zeroth harmonic mode), 
due to the asymmetric form of the nonlinear interaction potential. The existence of asymmetric 
envelope localized modes is predicted. The types and characteristics of these coherent nonlinear 
structures are discussed. 
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I. INTRODUCTION 



A remarkable new feature associated with the physics of dust-contaminated plasmas P, 01 is the existence of 
'"^ , strongly-coupled charged matter configurations, which even lead to the formation of dusty plasma (DP) Wigner-type 
^ layers (crystals) when the inter-grain potential energy far exceeds the average kinetic energy of the dust particles. 
O ] Dust crystals display a variety of new phenomena, such as phase transitions (crystallization, melting) and lattice 
excitations; a link is thus established between plasma physics and solid state physics 

Dust lattices, which are typically formed in the sheath region in discharge experiments and remain suspended above 
the negative electrode due to a balance between the external electric and gravity forces are known to support 
^ ■ harmonic oscillations (acoustic modes) in both longitudinal and transverse-shear (horizontal-plane) directions, as well 
', as optical-mode-like oscillations in the vertical (off-plane) direction [E IE S B ^| ■ 

Dust-lattice waves (DLWs) are reminiscent of waves ('phonons') propagating in atomic chains, which are long 
, known to be dominated by interesting nonlinear phenomena (localized modes, instabilities) , due to the intrinsic 
■ nonlinearities of inter-atomic interaction mechanisms and/or on-site substrate potentials [III [l3 . IT3. IT^ . The present 
' study is devoted to the study of one such phenomenon, namely the nonlinear amplitude modulation of weakly nonlinear 
oscillations with respect to longitudinal lattice oscillations. Even though this mechanism, which is associated with 
' harmonic generation due to self-interaction of the carrier waves, is today known to be of relevance in phenomena 
, as diverse as the modulational instability and energy localization in solids, charge and information transport in 
^ ■ biomolecules and DNA strands, and coherent signal transmission in electric lines and in optical pulse propagation 
nJLt ! 0. 0) ' it only recently been considered in dusty plasmas 1^ with regard to the propagation of 

^ • dust-acoustic and dust-ion acoustic waves [ll Q] in a weakly-coupled (gas-like) DP. As far as the nonlinear modulation 
Q of dust-lattice oscillations is concerned, our knowledge is still at a very early stage, let alone some first attempts to 
Q model the modulation of longitudinal and transverse |^ DLWs. In addition to a theoretical investigation of 
the occurrence of modulational instability, these studies have provided a first prediction of the possible existence of 
localized envelope excitations in a dust-crystal. It should be stressed that, as long known and extensively studied 
in solid state physics, the physics of these localized modes (basically associated with the generation of harmonics, 
M which was recently observed experimentally; see in Ref. |22|) is essentially different from that of compressive solitons, 
described by Korteweg-deVries (KdV)- type mod els, first introduced by Melands0 and later considered in some 
experiments [Ull^ and in theoretical papers [2 H25I l2l l27 | . 



In this paper, we consider the amplitude modulation of longitudinal waves propagating in DP lattices. We will 
derive an equation governing the evolution of the dust lattice wave envelope and will then examine its stability with 
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respect to external perturbations. Explicit conditions for the formation of different types of envelope excitations 
will be presented, in terms of the physical parameters involved in DP crystals and, in particular, the existence of 
asymmetric envelope dust-lattice solitons is predicted. This study refers to the longitudinal grain motion in an 'ideal' 
one-dimensional DP crystal, i.e. a single, unidimensional, infinite-sized, dust-layer of identical (in size, charge and 
mass) dust grains situated at spatially periodic sites (at equilibrium). It should however be noted that real crystals 
formed in the laboratory are often two- or three-dimensional arrays 128| , which may host longitudinal waves such 
as the ones considered herein, but also transverse (shear) waves 3 ^ weW as dissipative nonlinear structures, as 
observed in recent experiments j22,i2j|. Two-dimensional (2D) lattices are known from Solid State Physics to possess 
an enriched profile of possible excitations, including e.g. oblique (with respect to the principal axes) waves, 2D 
breathers and vortex excitations |29j| . Therefore, this work is only a first step towards the investigation of nonlinear 
mechanisms in DP crystals: the mechanism put forward in this paper may be generalized to higher dimensionality; 
nevertheless, this is left for consideration in future work. 



II. THE MODEL 



Let us consider a layer of charged dust grains (mass M and charge q, both assumed constant for simplicity) of 
lattice constant vq. The Hamiltonian of such a chain is of the form 



where r„ is the position vector of the n— th grain; Unm{fnm) = Qipix) is a binary interaction potential function related 
to the electrostatic potential 4>{x) around the m— th grain, and = |r„ — is the distance between the n— th and 
m— th grains. We shall limit ourselves to considering the longitudinal (^ x) motion of the n—th dust grain, which 
obeys 

where E{x) = —d(t){x) / dx is the electric field; the usual ad hoc damping term is introduced in the left-hand-side of 
Eq. involving the damping rate due to dust-neutral collisions. A one-dimensional (ID) DP layer is considered 
here, but the generalization to a two-dimensional (2D) grid is straightforward. At a first step, we have omitted the 
external force term F^xt, often introduced to account for the initial laser excitation and/or the parabolic confinement 
which ensures horizontal lattice equilibrium in experiments p^ . 



A. Equation of motion 



Assuming small displacements from equilibrium, one may Taylor expand the interaction potential 0(r) around the 
equilibrium inter-grain distance Itq = \n — m|ro (between Z— th order neighbors, 1 — 1,2, ...), viz. 



1 d^'(h(r) 



I'l dr^' 
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where V denotes the degree of nonlinearity involved in its contribution: T = 1 is the linear interaction term, V — 2 
stands for the quadratic potential nonlinearity, and so forth. Obviously, 5xn = Xn — xi"' denotes the displacement of 
the n— th grain from equilibrium, which now follows 
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We have distinguished the hnear/nonhnear contributions of the first neighbors (lst/2nd lines) from the corresponding 
longer neighbor terms (3rd/4th lines, respectively). 

Adopting the standard continuum approximation, often employed in solid state physics 0, we may assume that 
only small displacement variations occur between neighboring sites, i.e. 



5xn±i = Sxn ± Iro 



du 
dx 
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where the displacement Sx{t) is now expressed by a continuous function u = u{x,t). One may now proceed by 
inserting this ansatz in the discrete equation of motion |(5J), and carefully evaluating the contribution of each term. 
The calculation, quite tedious yet perfectly straightforward (to be reported, in full detail, elsewhere 113), leads to a 
continuum equation of motion in the form 



il + VU~Vq Uxx 
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where the subscript x denotes the partial differentiation, viz. UxUxx ~ k^x)xl'^ and (ux) Uxx — iu^x^xji 



B. Characteristic physical quantities 

The definition of the 'sound speed', wq = "-^o.l associated with the longitudinal oscillation eigenfrequency wo,l, 



i=\ 



where iV is the degree of interacting site vicinity assumed: A'^ = 1 for first-neighbor interactions (FNI), iV = 2 for 
second-neighbor interactions (SNI), and so forth. The dispersion coefficient, involving the characteristic velocity 
is given by 
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We note that v\ = Vq/12 for = 1 (only). The quadratic force (or cubic potential) nonlincarity coefficient is 
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Finally, the cubic force (or quartic potential) nonlinearity coefficient is 
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'?o = i^r^^r'(M^'- (7) 
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Let us point out that the above definitions of the coefficients in Q are inspired by the Debye-Hiickel (Yukawa) 
potential of the form (jyoir) — {q/r) e^^/^'^ (whose odd/even derivatives are negative/positive), in which case they 
are defined in such a way that all of Uq, v\, pq and qo take positive values. In this case (and for A^ = 1 i.e. FNI), 
relations igl), (0), lO and yield 
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where we have introduced the lattice parameter k = r^/Xi 
spontaneously formed in current dusty plasma dischar ges |^ 
coincide with the ones in previous studies for FNI |^ |23| ■ 



(15), (16), (26) in Ref. [|], 
the notation of Amin et al. 



>; in fact, k is roughly between 1 and 1.5 in DP chains 
IHEI^' These relations, which are depicted in Fig. ^ 
Specifically, compared to the notation of Melands0 [see 



and po above correspond to Wq, Wq/12 and j{a)a'^/M therein. Also, compared to 



[see (l)-(4) in Ref. 
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However, the coefficient qq, which can easily be seen to be practically 
to be omitted ~ is neglected in those studies. 

One should retain, for later reference, the form of the discrete dispersion relation 
(discrete) equation of motion Q 
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One may readily verify that the standard ID acoustic wave dispersion relation lu kvo is obtained in the small k 
(long wavelength) limit: check by setting sin(/ kro/2) « lkro/2 (and recalling the general definition of vq above). For 
FNI and Debye-Hiickel interactions, one obtains 
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[cf. for iV = 1]. We note that the LDLW frequency and the sound speed will be of the order of: Cjq — 
[2g2/(MA|))]^/2 (= clIo\/2 = LUmax) and vq = [2q'^ / {M Xd)]^^"^ = woA|,, respectively; these characteristic quanti- 
ties scale as ujo ^ r^^l"^ a^'^l'^ U and (I^o ~ r^^/^ a^^/^ U , respectively, for a given dust grain radius r, an inter-grain 
separation a and the dust surface potential U = q/r (following Ref. J5s|)- 

Even though all the formulae above are provided in their general form, in the following we shall assume that = 1 
(FNI) and will also set ^ « at a first step, for the sake of analytical tractability. These assumptions will be relaxed 
in a forthcoming work, also accounting for higher-order nonlinearities neglected here, which may be evaluated in a 
systematic manner 30] . Our aim here is to emphasize on a physical mechanism, rather than putting forward a tedious 
and exhaustive mathematical model. 



III. MULTIPLE SCALE EXPANSION - DERIVATION OF A NONLINEAR SCHRODINGER 

EQUATION 

According to the standard reductive perturbation method [Til iBlj , we shall consider a small displacement of the 
form: u ^ + eui -\- e^U2 -t- ... , where e ^ 1 is a small parameter, and solutions u„ at each order are assumed 

to be a sum of m— th order harmonics, viz. u„ = X]m=o ^ exp[i(fca:: — ujt)] (the reality condition uf?2t = '"m'' is 
understood). Time and space scales are accordingly expanded as 



and 



implying that 



and 



d/dt ^ d/dTa + e d/dTi + e^d/dT2 + ... = So + e 9i + e'Oa + .. 

d/dx dIdXo + e d/dXi + t^djdX^ + ... = Vq + e Vi + V2 4- 

d'^ldt^ ^ dl + 2edodi + e^{df + 28082) + ... , 
dyOx"^ Wl + 2e VoVi + e^Wj + 2V0V2) -f ... , 

d^/dx"^ -f 4e VjjVi +2e'^Wl (SV^ -f- 2V0V2) + ... . 



We now proceed by substituting all the above series in (j2Jl and isolating terms arising in the equation of motion at 
each order in e". By solving for u„, then substituting in the following order e"+i and so forth, we obtain the m— th 
harmonic amplitudes ui'"'' in each order, along with a compatibility condition up to any given order. 
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The equation obtained in order is 

id^-vlVl~vlrlvl)ui = Loui =0. (13) 

We may assume that ui — uj^' exp[i{kx ~ ut)] + c.c. = exp{i6) + c.c., where u, k — 2n/X and A denote the carrier wave 
frequency, wavenumber and wavelength, respectively; c.c. stands for the complex conjugate of the preceding quantity 
everywhere. 

The dispersion relation obtained from H13() is of the form 

2 , 2 2 ; 4 2 2 ; 2 2 '^1 ; 2 2\ 

uj = k Vq ~ k v^Tq — k Vq yl — —k Vq j , (14) 

predicting an acoustic behaviour (viz. uj « kvo) for very low k, and dispersion (weighted by w^J'o) higher values. 
Notice that the reality condition a; > 0, implying k < (vo/vi)rQ^ = kcr,i (= V^^r^^ for FNI), is in principle 
satisfied, covered by the continuum hypothesis k <C (7r/ro) = kcr,o- Specifically, one may numerically verify the 
agreement between the discrete dispersion relation in (|12|) (for v — 0) and its continuum approximation (|14|l - see 
Fig. 121 for instance, aiming at an accuracy interval of 99 % (90 %), i.e. allowing for a relative error of 1 % (10 %), 
one should restrict the confidence interval of this study to values of fcrg below 0.3 (1.0). 

Let us evaluate the action of the linear operator Lq, defined above, on higher harmonics of the phase 9 

Lo e"^ = [(-^nw)2 - v^^ink)^ - r^iinkf] e"^ = ... = -n\n^ - l)vl k* e"^ = D„ e™^ , (15) 

where we made use of the dispersion relation (|14|) . We observe that Dq ^ Di = 0. 
To order ~ e^, we have 

Lo U2 = -2{do di - wgVoVi - 2vl rg Vj^Vi) ui - po {VqUi) (V^ui) . (16) 

One now has to impose the condition of suppression of secular terms (i.e. terms ^ e'^ in the right-hand side, which 
may resonate with the null space of the operator Lq), which takes the form 



dTi " dXi 



= 0, (17) 



or du'^i^ /dC = 0, i.e. u'^i^ — u^i\^,t), where C — e{x — Vgt) and r = e^t, implying that the slowly-varying wave- 
front (envelope) travels at the group velocity Vg = uj' (k) — (wq — 2vir'^k'^)k / uj . We notice that Vg becomes zero at 
k = uo/(ui\/2) r-Q ^ = fccr, i/^/2 \/6rQ^ for FNI), beyond which it becomes negative: the wave envelope and the 
carrier wave then propagate in opposite directions (backward wave). Nevertheless, for the sake of rigor, note that this 
comment rather becomes obsolete once one properly takes into account the influence of the lattice discreteness on the 
form of the dispersion relation uj = uj{k); see Fig. |21 

(2) 

The solution U2 obtained at this order involves a second harmonic contribution U2 exp{2i6) = 

{ipok'^ / D2){u^i''Y exp(2i0) = —i[pQ / {12v\rQk)\ (u^^^)^ exp(2i0), as well as a zeroth (and possibly a first) harmonic(s), 
which may be determined by higher-order relations. 

One may now proceed to the next (third) order in e. The condition for annihilation of secular terms now takes the 
form 

The dispersion coefficient P is 

r, 1 ^ t "2 2 , a 2 2i2\ ^1^0^^ ( i2 ^Up \ 

^=2" ('')=-^(^9--o+Qvirok)^^^[k (19) 



while the nonlinearity coefficient Qq is 
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For qo = and Ug^'' = 0, these expressions agree exactly with Eqs. (17) and (18) in Ref. '2^. However, the zeroth- 
harmonic correction Wq^^ cannot be set equal to zero; indeed, collecting the constant (zeroth-harmonic) terms arising 
in this order, one obtains 

(21) 

which may immediately be integrated, yielding 

^^0^^ Pofc^ L.(l)|2 , _ 73^i.M.,(l)|2 



dc 



C = R{k)\uY>\' + C, (22) 



where C is a constant, to be determined by the boundary conditions (e.g. C = if m^^^ = at a; = ±(X)). Note that 
the resonance condition Vg — vq amounts to k'^ — 3tig/(4w^rQ) = 0, which is not satisfied in the parameter range of the 
continuum approximation; on the contrary, Vg < vo everywhere, implying that R{k) takes positive/negative values for 

positive/negative pq: for Debye interactions, po > and R > 0. The non-vanishing constant term Uq^^ is purely due 
to cubic interaction potential nonlinearity (and disappears for po = 0). This is known in solid state physics jTH . 

Eqs. H18|l and H22(l form a coupled set of evolution equations, governing the dynamics at first order in e. They may 
readily be combined into the form of a Nonlinear Schrodinger Equation (NLSE) 

^^ + P^ + Q\A\^A^O, (23) 

which, once solved for the amplitude of the first harmonic A — u^i \ immediately provides the solution for u'^'^ 
upon substitution into H22() . A linear term ^ C A has been neglected, since it is canceled by a trivial phase shift 
transformation. The final form of the nonlinearity coefficient Q now reads 

Q = -Tr[iok^ + 7rh) ~ — 2- 24 

In the infinite wavelength limit (i.e. near fc w 0), we have 

p^_^]±lk~^-^k^ + 0{k^) 



and 



Let us stress that the extra term, which makes the difference between Qo and Q, changes completely the behaviour 
of the nonlinearity coefficient near fc w 0. For example, we find that 

Note, once more, that the results in this Section are general and apply to a wide class of physical problems. In 
any given problem, involving an arbitrary (analytical, possibly long-range) interaction law, Eqs. (|18l) . (|22|1 and H23() 
provide the evolution of the first-order correction ui to the particle displacement u, if the latter is governed by an 
equation of motion of the form of (j^J, as expressed in the continuum limit by Eq. (PJ. To first order in e, the solution 
then reads 

u{x, t) « e [liW r) + it« (C, r) + c.c] + 0{e^) , 

where u^f^ is given by \22\ (it vanishes if, and only if, po = 0) and u'-p corresponds to (the amplitude of) a small 
translation and an oscillation, respectively, both traveling at the group velocity Vg. We emphasize, as a general result, 
that the existence of a non-zero cubic potential nonlinearity term po implies the existence of a finite translation 
(zeroth-harmonic) correction Uq^' to first order in e, while the sign of pq affects the form of u^^^ (i.e. the lattice 
compression or rarefaction); furthermore, it also determines the (sign of the) second-harmonic correction, as we saw 
above. Nevertheless, the sign of po does not affect the form of the coefficients in (I23|l . and thus neither the plane wave 
stability profile. 
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IV. STABILITY ANALYSIS 



The standard stability analysis [iSl 133 . l33l| consists in linearizing around the (Stokes) plane wave solution of l|23|l . 
A = Ae*'^!"^!^'^ + c.c. (notice the amplitude dependence of the frequency), by setting A ^ Aq + eAi and taking the 



perturbation Ai to be of the form Ai 



c.c.^ where k and tj denote the perturbation's wavenumber 



and frequency, to be distinguished from the carrier wave's homologous quantities, k and lo. Hence from (|23f) one 
readily obtains the dispersion relation di'^ — k"^ ( k^ — 2^\Ao\'^^ . The wave will thus be stable if the product PQ 

'2^\Ao\, 



is negative. However, for positive PQ > 0, instability sets in for wavenumbers below a critical value kcr — y ^pi 
i.e. for wavelengths above a threshold: Ac- = 2'K/kcr, with an increment cr = \Imuj{k)\ which attains a maximum 



value at A; = kcr/V^, viz. 



— llmuj] 



k=k^ 



= \Q\\M 



It turns out that the instability condition depends only on the sign of the product PQ, which may be studied 
numerically by using the expressions derived above. 

Let us go back to the expressions p9ll and (|24|l for the coefficients P and Q. First, we notice that P is always 
negative, given the condition k < kcr,i = (vo/vi) r^^ imposed above. On the other hand, the nonlinearity coefficient 
Q may take either negative or positive values, depending on the physical parameters vq, ui, po and qo and, of course, 
the wavenumber k. To see this, one may cast (|24f) in the form 



Q 



"2uj 



qo k 



3 - 2vlrlk^/vl 



(25) 



and then examine its sign in terms of k. This can be done either by a straightforward (yet rather lengthy) analytical 
calculation (omitted here) or numerically (see below). Notice, as a general remark, the pole created at the resonance 
point Vg = wq, for k = VQ^/i/{2vi)r'^^ (= ir^^ for FNI). In conclusion, the carrier wave will be modulationally 
unstable for k values corresponding to Q < [PQ > 0), and stable for Q > {PQ < 0). These results are valid for 
any given form of the electrostatic interaction potential (p{r) and neighboring site vicinity N, which (both) in fact 
enter the physical parameters via the expressions - 10). 

Let us now focus on the dust-lattice problem. At a first step, one may explicitly assume first neighbor Debye- 
Hiickel type interactions, by substituting with expressions © - H1U|) into the above relations for P, Q, and then study 
the behaviour of all the relevant physical quantities as a function of, say, the lattice parameter k. Notice that the 
stability profile deduced above is dramatically modified, in comparison with the limit of vanishing cubic interaction 
nonlinearity, i.e. qo — 0, considered in Ref. |20|. where the instability of plane waves was directly prescribed (in 

addition to Mq^^ =0)- A simple numerical investigation shows that Q is always positive (prescribing stability, since 
P < 0) at low values of the wavenumber k; see Figures 13 0] This is true everywhere, for k below the resonance 
threshold kro « 3, except for k of the order of « 0.75 r^^ or lower, where a small instability region is encountered; see 
Fig. 0] [notice the asymptotic curve at kro = i;oV3/(2ui) = 3, as expected from H25|l. beyond which the continuum 
approximation certainly fails]. Long wavelength LDLWs are therefore expected to be stable, while shorter wavelength 
LDLWs may be unstable. For rigor, it should be admitted that the instability (carrier wavenumber) threshold kcr 
is quite high, and rather beyond the validity of this continuum model (recall the assumption made above, that 
k <C kcr,o = 7r/r0). Therefore, while small k stability may be taken for granted, we only possess an indication for (i.e. 
not a rigorous proof of) instability for longer k. We draw the conclusion that LDL waves are most likely to remain 
stable with respect to external perturbations (even though this statement may a priori not be valid for a different 
interaction potential). The present results are in agreement with previous prediction in atomic chains [TllIT^ . 



V. NONLINEAR EXCITATIONS 

Equation H23|l is an integrable nonlinear partial differential equation [33 . l34j , which is known to possess different 
types of localized, constant profile, travelling wave solutions (solitons). Several types of such localized modes are 
presented in Refs. |0, 0, ) (see Ref. for a review) so we will only briefiy outline the analytical form of 
those we are interested in, and discuss their relevance to our problem. 

A solution of Eq. (|23|l may be sought in the form A{(,t) = piCi''') e*®^^''^^ where p and a are real variables to be 
determined. Different types of solution are thus obtained, depending on the sign of the product PQ. Once wj^^ — A 
is thus determined, the zeroth-harmonic first-order correction Uq^"* is obtained from (|22|l . 
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For PQ > we find the bright envelope soliton 

A{x,t) = pq sech[{x — v^t) / Lf, + Xe] a{x — Vgt) — 4,{a^ — rf) P t] (26) 

{sechx = 1/ cosh a;), which represents a localized envelope pulse travelling at a speed — Vg — AaP. The pulse 
width Le — l/{2f]) is related to the (constant) maximum amplitude po = '^viP/Q)^^^ ^is Po = (P/QY^^ /Le- See that 
both of the small (arbitrary) parameters a and rj bear dimensions of inverse length. Substituting H2t)|) into H22|) , we 

obtain the zeroth-harmonic term u'^^ = i? J dC = 2r]PR{k)/Q tanh[(a; — Vet)/Le + Xe] + c (we will assume 

that c = 0). Combining the preceding formulae, we obtain the asymmetric bright-envelope modulated solution for 
the (total) dust-grain displacement (to order ~ e) 

Ml = ^0 sech[{x — Vf.t)/Le + Ce] cos[(a:: — Vct)/ Lq + Cq] + Bq tanh[(a; — v^t)/ + Cg] (27) 

where Aq — 2po = 477(P/Q)^/^, Bq = R{k)p^Le — 2riR{k)P/Q and the shifted (modulated) carrier wavenumber 
kc = Lq^ and frequency ujc — vq/Lq now depend on: Lq — \/(k — 2a) and Vc — [uj — 2aVg + 4P (a^ — ri^)]Lo; 
Ce and Ce are arbitrary real constants. Equation (|27|l represents a superposition of a localized modulated harmonic 
oscillation and a kink-like excitation, which is characterized by different constant asymptotic limits (at x — ±oo). 
Since R(k) > in a Debye crystal, this solution will therefore represent a localized compression co-propagating with 
a strongly modulated envelope oscillation; see Fig. E 
For PQ < 0, we have the grey envelope soliton [3g 

^(C, r) - pi {1 - sech^iiC - 2ar)/Li]}i/2 exp[z aiC, t)] , (28) 

where 

a{C, t) = sin"i{atanh(C/Li) [1 - sech'^{C/Li)]-^/'^ - 2aC + {Q - iPa'^^} 

which represents a localized region of negative wave density (a void), with finite amplitude (1 — a)pi at ^ — 0; 
< a < 1) [one may choose fl = iPrj^ and pi — 2 77 (jP/Qj)^/^, for analogy with l|26|l ]. Again, the pulse width 
Li = {\P/Q\y/^/{api) is inversely proportional to the amplitude pi. Going back to the original variables, we obtain 
the dust grain harmonic oscillation 

u[^^ expie + c.c. = Ai{l- c? sech^\{x ~ Vet)/Li + Ce]}^/^ cos[(a; - Vct)/Lo + ai{x,t)] , (29) 

which has a maximum amplitude Ai = 2 pi and a phase shift 

ai{x, t) — sin^"'"[atanh(a:: — Vet)/Li] [1 — sech^lx — v^t)/ Li]~^^'^ + (Tq , 

where Ug, Vc and Lq are just as defined previously. Notice the (dimensionless) parameter a, which regulates the depth 
of the excitation. For a = 1, one obtains the dark envelope soliton 

u'^i ^ e'xjpid + c.c. = Ai tanh[(a; — Vct)/Li + Cg] cos[(a; — Vct)/Lo + cq] , (30) 

which describes a localized density hole, characterized by a vanishing amplitude at C = 0; see Fig. The dark- 
or grey-soliton-modulated wavepackets in the previous paragraph, are superposed on a constant (zeroth-harmonic) 
displacement, given by as 

Wq^-* = —Bi ta.nh[{x — Vet) / Li + Ce] , (31) 

where Bi = R{k) p^a^ Li = aRik) (|P/Q|)^^^ pi; the constant C in (|22(l was chosen as C = —R{k)p\ for constant 
asymptotic values at infinity. Combining the latter formulae, we obtain the asymmetric dark/grey modulated envelope 
solution for the dust-grain total displacement (to order '-^ e) 

m = Ai{l- sec/i^[(C - 2ar)/ii]}^/^ cos[(a; - Vet)/Lo + a + cq] - Bi tanh[(a; - Vet)/Le + Ce] (32) 
or, for a ~ I, 

ui — Ai tanh[(x — Vet)/Li + Ce] cos[(x — Vct)/ Lq + cq] — Bi tanh[(a; — Vet)/Le + Ce] (33) 

where all quantities were defined previously. For i? > 0, as in the Debye DP crystal case, the dark/grey type excitations 
correspond to a localized density dip accompanied by a rarefaction propagating in the lattice; see Fig. |3 
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A quantity of importance comes out to be the ratio ^ = P/Q: its sign dictates the type of envelope soliton which 
should, in principle, propagate in the lattice, while (the square root of) its absolute value represents a measure of the 
soliton amplitude for a given width L, viz. Aq = \/^/L. In figure [7| we have depicted ^ versus the (normalized) 
wavenumber kro and the lattice parameter k. We notice an increase of the absolute value of the (negative, mostly) 
ratio with decreasing k and with increasing k. One draws the conclusion that for a given excitation width L (assumed 
L 3> To, according to our approximation), shorter wavelength (i.e. higher wavenumber) carrier waves will form 
dark-type excitations (pulses) with an increased maximum displacement Aq] this picture, which is true below, say, 
k « 2.25 Vq^, is inverted above this value (nevertheless, there may be some doubt above the validity of results in 
that region, as discussed above). Also, higher values of k, i.e. higher lattice constant rg values for a given Debye 
length Xo, seem to favor narrower excitations (as compared to tq). Notice, in Fig. {7^, the curly-shaped curve at the 
bottom, indicating the poles of the ratio P/Q, i.e. where Q = (i^ > only below this curve), as well as the resonance 
asymptotic line at kro = 3. See that positive values of ^ are confined in a small, rather short-wavelength and small k 
region; see Fig. [7^. As a final comment, notice that values deduced from Fig. [7|are in agreement with the continuum 
approximation. For instance, as a crude estimation, taking L = 5ro and k « 1, one obtains ^ ~ 0.025rQ, implying 
^0 = V^/L « O.OSro, i.e. Aq <^ ro as predicted. 

Concluding, the instability/stability regions depicted e.g. in Figs. O ^ in fact also limit the (fc, k) parameter 
pair values where bright (dark/grey) solutions i.e. density pulses (holes) may exist. Furthermore, the envelope 
characteristics will depend on the carrier wave dispersion via P and Q; for instance, regions with lower values of P 
(or higher values of Q) will support narrower excitations. Admittedly, however, dust lattices will rather not favor 
bright-type structure formation by unstable long-wavelength LDL waves. On the other hand, such waves will mostly 
be stable to external perturbations, and may locally form density dips (envelope holes). 

VI. CONCLUSIONS 

In this paper, we have investigated the amplitude modulation of longitudinal dust-lattice waves propagating in 
a one-dimensional dust-crystal. By using a continuum approximation, we have derived a Nonlinear Schrodinger 
Equation governing the evolution of the LDL wave envelope. We have found that long wavelength carrier waves 
will generally be stable, possibly evolving towards the formation of voids, i.e. asymmetric dark/grey-type coherent 
envelope structures, while longer wavelengths may allow for the formation of envelope pulses (asymmetric bright- 
type envelope excitations). Furthermore, we observe that the total grain displacement will bear a (weak) constant 
displacement (zeroth mode), due to the intrinsic third-order nonlinearity of the Debye potential. It may be admitted 
that the former (dark) excitations may rather not be relevant to this idealized, infinite chain model, since they 
correspond to an infinite amount of energy stored in the lattice. Nevertheless, they may be taken as a hint towards 
excitations sustained in finite-length, highly discrete, real lattices. 

The present study aims in making a first step towards the elucidation of the nonlinear harmonic generation mech- 
anisms dominating weakly-nonlinear longitudinal wave packets in dust-crystals. The formation of LDLW-related 
structures has thus been predicted, and will hopefully be confirmed by appropriately designed experiments or by 
numerical simulations. In fact, through the present investigation, a link has been established towards the formalism 
of discreteness-related intrinsic localized modes (discrete breathers) widely studied recently in solid state systems 
[stLIs^ . This work should, in principle, be extended by including realistic effects associated with crystal asymmetries, 
defects, dust charging, ion-drag and multiple dust-layer coupling 39]. In particular, one should consider a similar 
effect in more realistic, e.g. 2D, geometries, as the ones observed in recent experiments with laser-excited dust lattices 
|22ll24j . Concluding, this work was devoted to a nonlinear mechanism for harmonic generation and coherent structure 
formation. Although here presented in a ID geometry, for simplicity, this qualitative mechanism may be generalized 
to 2D lattices; this possibility will be considered in future work. 
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Figure captions 

Figure 1. 

(a) The characteristic physical quantities defined in - poil (for Debye-type first neighbour interactions) are 
depicted against the lattice parameter k: the eigenfrequency (square) a;£ g, normalized by Q'^/MX^ (solid curve — ); 

the sound speed (square) Vq, normahzed by Q'^/MXd (long-dashed curve ); the cubic potential nonlinearity 

coefficient po, normahzed by Q^/MXd (dot curve • • •); the quartic potential nonlinearity coefficient go, normalized by 

/{MXdY/^; short-dashed curve ); (b) close-up near k « 1. 

Figure 2. 

(a) The linear oscillation frequency uj (normalized byr — {2Q^ /MX^Y^'^) is depicted against the normalized 
wavenumber kr^; the three curves, from bottom to top, depict: the exact relation Hll|l ( — ), the continuum approxi- 
mation H14|l ( ), and the (tangent) acoustic limit lu = kvQ ( ); (b) close-up near the origin. 

Figure 3. 

The NLS coefficients P, Q and their ratio ^ ~ P/Q'- discrete model [see (|11() ] (dashed curve) vs. continuum model 
(straight curve), are depicted vs. the (normalized) wavenumber fcrg, for k = 1.25. (a) P, cf. 119|l (normalized by 
(Q2A£,/M)1/2); (b) Q, cf. ^ (normalized by [QV(MA^)]1/2); (c) C = P/Q (normalized by X%- divide by a factor 

— 1.25'* K, 2.44 to obtain the value in rg). 

Figure 4. 

Similar to the previous figure, for k = 0.25. The continuum model seems to fail completely above k = \: notice the 
discrete (correct) curve, which is hardly visible (dashed curve). 

Figure 5. 

(a) Bright (asymmetric) envelope solution {PQ > 0): the first- and zeroth- harmonic components of the particle 

displacement, u^^-* and u'^\ are depicted, for t = 0; dash ( ) and solid ( — ) curves, respectively, (b) The total 

displacement ui (~ e) is depicted (arbitrary parameter values). 

Figure 6. 

Asymmetric envelope solutions for PQ < 0: the total displacement is depicted for t — (arbitrary parameter 
values): (a) grey-type envelope solution (0 < \a\ < 1); (a) dark-type envelope solution (|a| = 1). 

Figure 7. 

Contours of the reduced ratio ^ — P/Q (expressed in units Tq) are depicted vs. the (normalized) wavenumber kro 
and the lattice parameter k. (a) The ^ values are, from top left and downwards: —0.002, —0.005, —0.01, —0.02, —0.03. 
The absolute value of the (negative) ratio increases with decreasing k and with increasing k. Notice the curly-shaped 
curve at the bottom, indicating the poles of the ratio P/Q, i.e. where Q = (^ > below this curve), as well as the 
resonance asymptotic line at kro — 3. (b) Close up for small k and positive ^ (wave instability region). The ^ values 
are, from bottom to top: 0.35, 0.4, 0.5, 1, 2 (and -l-oo, again, i.e. Q = 0, on top). 
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